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r '~j ' Abstract. The GRB010222 optical and near-infrared (NIR) afterglow was monitored at the TNG and other Italian telescopes 

K> ' starting ~1 day after the high-energy prompt event. The BVR light curves, which are the best sampled, are continuously 

^ . steepening and can be described by two power laws, f{t) cx t~", of indices ai ~ 0.7 and Q2 ~ 1-3 before and after a 

■ - - ' break occurring at about 0.5 days after the GRB start time, respectively. This model accounts well also for the flux in the 

U, I and J bands, which are less well monitored. The temporal break appears to be achromatic. The two Tf-band points are 
not consistent with the above behaviour, and rather suggest a constant trend. A low-resolution optical spectrum has also been 
taken with TNG. In the optical spectrum we found three absorption systems at different redshifts (0.927, 1.155 and 1.475), 
the highest of which represents a lower limit to, and probably coincides with, the redshift of the GRB. The broad-band optical 
spectral energy distributions do not appear to vary with time, consistently with the achromatic behaviour of the light curves. We 
compare our measurements with different afterglow evolution scenarios and we find that they favor a transition from relativistic 
to non-relativistic conditions in the shock propagation. 
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Nazionale Galileo (TNG), operated on the island of La Palma by 
the Centro Galileo Galilei of the CNAA (Consorzio Nazionale per 
I'Astronomia e 1' Astrofisica) at the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Canarias; the 
Asiago Astronomical Observatory, Italy; the Bologna Astronomical 
Observatory in Loiano, Italy; the Campo Imperatore Astronomical 
Observatory, Italy, and; TIRGO infrared observatory, Switzerland 



1. Introduction 

The light curves of Optical Transients (OTs) associated with 
GRBs are generally described by single power laws f{t) oc 
i^" with indices a ~ 1.1 — 2. However, in a number of well 
monitored cases a change in the light decay rate at about 0.5-1 
days after the GRB was detected with a transition to a steeper 
power law behaviour. The power law indices before and after 
the temporal break are typically in the range 0.7-1 and 1.7-2.4, 
respectively (GRB990123: Fruchter et al. 1999, Castro-Tirado 
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et al. 1999, Kulkarni et al. 1999; GRB990510: Stanek et al. 
1999, Harrison et al. 1999, Israel et al. 1999; GRB990705: 
Masetti et al. 2000a; GRB000926: Fynbo et al. 2001a). This 
behaviour can be caused by the deceleration of a relativistic jet 
in a homogeneous medium (Sari et al. 1999, Rhoads 1999), 
or by expansion in a wind (ChevaUer & Li 1999, 2000), or 
by the transition between relativistic and Newtonian conditions 
in the plasma expansion (Dai & Lu 1999). Fits to multiwave- 
length data of individual afterglows with the above models have 
been proved to be satisfactory (Dai & Lu 1999, Panaitescu & 
Kumar 2000). The multiwavelength optical and near-infrared 
(NIR) spectra of OTs are generally well accounted for by syn- 
chrotron radiation in a relativistically expanding shock, and are 
described by power laws of different characteristic indices, de- 
pending on the location of injection and cooling breaks, and 
temporal evolution thereof (Sari et al. 1998). Both the fading 
rate and the spectral slopes are determined by the shape of the 
electron distribution, in different ways, according to the geom- 
etry of the emitting region. Departures of the optical-NIR spec- 
tral slopes from those expected based on the flux temporal de- 
crease are often interpreted as due to absorption of the OT light 
within the ORB host galaxy (e.g. Palazzi et al. 1998, Dal Fiume 
et al. 2000, Klose et al. 2000, Price et al. 2001a). 

GRBO 10222 was simultaneously detected by the Gamma- 
Ray Burst Monitor (GRBM) and the Wide Field Camera 
(WFC) unit 1 onboard BeppoSAX on 2001 Feb 22.30799 UT 
(Piro 2001a,b) as one of the most intense GRBs observed by 
both the GRBM and WFC. A quick repointing of the satel- 
lite allowed the detection of the X-ray afterglow at a position 
consistent with that of the prompt event (Gandolfi 2001). A 
detailed analysis of the high-energy data is presented by in 't 
Zand et al. (2001). 

Many observers started a ground-based campaign in order 
to search for the GRB afterglow at lower wavelengths. A bright 
transient source {R ^ 18.4) was independently detected in the 
optical by Henden (2001a,b; see also Henden & Vrba 2001) 
and McDowell et al. (2001) about 4 hours after the GRB. The 
object, which is not present in the DSS-II sky survey, is at co- 
ordinates (J2000) a = 14'' 52'" 12!55, 5 = +43° 01' 06'.'2 with 
an error of 0'.'2 along both directions (Henden & Vrba 2001), 
well inside the error box of both high-energy prompt event and 
X-ray afterglow as detected by the BeppoSAX Narrow Field 
Instruments. The OT is among the brightest ever observed as- 
sociated with GRBs. 

The fading behaviour of the object, first reported by Henden 
& Vrba (2001) and later confirmed by Stanek et al. (2001a,b), 
left little doubt on its afterglow nature. Assuming for the R- 
band flux a power law decay, typical of optical afterglows. Price 
et al. (2001b) and Fynbo et al. (2001a) measured a temporal 
slope a = 0.89±0.09 and a = 0.86±0.01, respectively. Masetti 
et al. (2001) reported that a steepening of the decay slope to 
a ~ 1.3 occurred at about 1 day after the GRB start time. 

Garnavich et al. (2001) acquired an optical spectrum of the 
OT soon after its detection. From the presence of Fe ii and 
Mg II absorption features, they measured a redshift z = 1.477 
for the afterglow. Jha et al. (2001a,b) also noticed the presence 
in the same spectrum of a second absorption system, located 



at z = 1.157. A third absorption system at z = 0.928 in the 
optical spectrum of the afterglow was reported by Bloom et al. 
(2001). The farthest of the three systems appeared to be formed 
by substructures at sUgtly different redshifts, consistent with 
internal gas motions in galaxies (Castro et al. 2001). 

A bright counterpart was also detected at NIR (Di Paola et 
al. 2001), submilhmetric (Fich et al. 2001) and radio (Berger & 
Frail 2001) wavelengths. 

The good sampling and the optical brightness of the 
GRBO 10222 afterglow have allowed a detailed study of its evo- 
lution at least up to about 40 days after the explosion. In this 
paper we present the results of the optical and NIR monitor- 
ing of the optical/NIR transient associated with GRB010222 
conducted at various Italian telescopes and started about 1 day 
after the GRB. Section 2 illustrates the observations and the 
analysis of the photometric and spectroscopic data; the results 
are presented in Section 3 and discussed in Section 4. 

2. Observations and data reduction 

We observed tiie GRBO 10222 field at the 3.58-meter Telescopio 
Nazionale Galileo (TNG) in the Canary Islands (Spain), at 
the 1.8-meter "Copernicus" telescope at Cima Ekar of the 
"Leonida Rosino" Astronomical Observatory of Asiago (Italy), 
and at the University of Bologna 1.52-meter "G.D. Cassini" 
telescope in Loiano (Italy). 

TNG was equipped with the spectrophotometer DOLoReS 
and a 2048 x 2048 pixels Loral CCD which allows covering a 
9.'5 x9f5 field in imaging mode with a scale of O'.'275/pix; the 
"Copernicus" telescope was carrying the AFOSC instrument 
whose 1100x1100 pixels SITE CCD has a field of view of 
8.'5 x8f5 with a scale of O'.'47/pix; the "Cassini" telescope was 
mounting the spectrophotometer BFOSC with a 1340x1340 
pixels LEV CCD, a field of view of 12^5 x 12^5 and a scale of 
O'.'58/pix. 

NIR imaging in J and K bands was acquired at the 
AZT-24 1 . 1 -meter telescope at Campo Imperatore (Italy) with 
SWIRCAM, a 256x256 pixels infrared camera with a field 
of view of A'AxA'A and a pixel scale of l'.'03/pix, and in 
J and Kg at the 1.5 -meter Gornergrat infrared telescope 
TIRGO (Switzerland) using the 256x256 pixels infrared cam- 
era ARNICA which has a 4' x4' field of view and a pixel scale 
of O'.'97/pix. The Kg filter is centered at 2. 12 jjm and has a full 
width at half maximum of 0.34 jim. For each NIR observation 
the total integration time was spUt into images of 30 s each, and 
the telescope was dithered in between. 

The complete log of our imaging observations is reported 
in Table 1. 

Two 30-min spectra, nominally covering the 3000-8000 
A optical range, were also acquired with TNG+DOLoReS be- 
tween Feb 23.227 and Feb 23.274 UT. The use of DOLoReS 
Grism #1 secured a dispersion of 2.4 A/pix. The sUt width was 
1'.'5. 

Optical images were debiased and flat-fielded with the stan- 
dard cleaning procedure. In some cases, frames taken on the 
same night in the same band were summed together in order 
to increase the signal-to-noise ratio. In Fig. 1 we report our 
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Table 1. Journal of the optical and NIR observations of the GRB010222 afterglow. Magnitude uncertainties are at la confidence 
level 



Exposure start 


Telescope 


Filter 


Total exposure 


Seeing 


Magnitude^ 


(Ui) 






time (s) 


(arcsec) 




2001 Feb 22.982 


AZT-24 


J 


2700 


4.0 


18.6±0.22 


23.056 


Asiago 


R 


1200 


3.8 


19.75±0.05 


23.063 


Asiago 


R 


600 


3.8 


19.79±0.04 


23.069 


Asiago 


I 


600 


3.8 


19.27±0.06 


23.082 


AZT-24 


K 


2400 


4.0 


17.2±0.32 


23.083 


Asiago 


I 


600 


3.8 


19.32±0.06 


23.173 


Loiano 


B 


2400 


4.5 


20.72±0.08=' 


23.199 


Loiano 


V 


1800 


4.0 


20.37±0.10^ 


23.204 


AZT-24 


J 


900 


4.0 


18.4±0.3^ 


23.211 


TNG 


R 


120 


1.1 


20.00±0.01 


23.219 


TNG 


R 


120 


1.1 


20.01±0.01 


23.280 


TNG 


U 


300 


1.1 


20.34±0.03 


23.284 


TNG 


U 


300 


1.1 


20.33±0.03 


24.097 


TIRGO 


J 


4680 


3.8 


19.21±0.35 


24.106 


Loiano 


V 


2400 


2.5 


21.23±0.06 


24.117 


TIRGO 




4680 


3.3 


17.49±0.25 


24.127 


Loiano 


R 


900 


2.0 


20.89±0.09 


24.236 


TNG 


R 


120 


1.0 


21.06±0.03 


24.239 


TNG 


R 


120 


1.0 


21.05±0.02 


24.241 


TNG 


I 


120 


1.0 


20.42±0.04 


24.244 


TNG 


I 


120 


1.0 


20.51 ±0.04 


24.247 


TNG 


V 


120 


1.0 


21.44±0.02 


24.249 


TNG 


V 


120 


1.0 


21.48±0.03 


24.252 


TNG 


B 


300 


1.0 


21.87±0.02 


24.257 


TNG 


B 


300 


1.0 


21.88±0.02 


24.262 


TNG 


U 


450 


1.0 


21.33±0.04 


24.268 


TNG 


u 


450 


1.0 


21.30±0.04 


25.253 


TNG 


R 


120 


0.9 


21.64±0.03 


25.261 


TNG 


V 


120 


0.9 


21.92±0.04 


25.263 


TNG 


V 


120 


0.9 


21.96±0.03 


25.267 


TNG 


B 


300 


0.9 


22.43±0.02 


25.271 


TNG 


B 


300 


0.9 


22.47±0.03 


25.275 


TNG 


I 


120 


0.9 


21.14±0.06 


25.277 


TNG 


I 


120 


0.9 


21.16±0.06 


27.139 


Asiago 


R 


900 


3.1 


22. 11 ±0.25 


27.264 


TNG 


R 


360 


1.2 


22.38±0.05 


27.270 


TNG 


V 


360 


1.1 


22.80±0.06 


27.278 


TNG 


B 


360 


1.1 


23.28±0.08 


Mai- 31.226 


TNG 


R 


3300 


1.2 


25.1±0.2 



Magnitudes of the GRB counterpart, not corrected for interstellar absorption 
^This measurement supersedes the value reported by Di Paola et al. (2001) 
^This measurement supersedes the value reported by Bartolini et al. (2001) 



second TNG i?-band image. We chose standard Point Spread 
Function (PSF) fitting as photometric method, and to this aim 
we used the DAOPHOT II image data analysis package PSF- 
fitting algorithm (Stetson 1987) running within MIDASQ A 
two-dimensional Gaussian profile with two free parameters 
(the half width at half maxima along x and y coordinates of 



' MIDAS (Munich Image Data Analysis System) 
is developed, distributed and maintained by ESO 
(European Southern Observatory) and is available at 

http : / / www .eso.org/projects/ esomidas 



each frame) was modeled on at least 5 unsaturated bright stars 
in each image. The errors associated with the measurements 
reported in Table 1 represent statistical uncertainties (at la) 
obtained with the standard PSF-fitting procedure. In only one 
case, i.e. the TNG observation of Mar 31, we used aperture 
photometry for the magnitude determination as the PSF-fitting 
procedure could not give reliable results due to the faintness of 
the transient. In this case we used an aperture diameter equal to 
the seeing FWHM of the summed image. 
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Table 2. List of the absorption lines identified in the TNG optical spectrum of GRB010222 afterglow. The number associated 
with each line refers to the identification shown in Fig. 4. The error on all line positions can conservatively be assumed to be 
±3 A, i.e. comparable with the pixel size (see text). EWs of lines in the ORB rest frame, i.e. divided by a factor (1+z), are also 
reported 



Line 


Observed 


Rest frame 


Element 


Redshift 




number 


wavelength (A) 


wavelength (A) 






(A) 


1 


3588 


1862.790 


Aim 


0.926±0.002 


2.1±0.2 


2 


3979 


2062.664 


Znll/Crll blend 


0.929±0.001 


2.5±0.3 


3 


4573 




Fell 


0.926±0.001 


0.4±0.2 


4 


4590 


2382.765 


Fell 


0.926=b0.001 


1.2±0.3 






ZOUU. L 1 J 


poll 


u.yzoztu.uui 


u. / =tU. J 


6 


5384 


2796.352 


Mgll 


0.925±0.001 


0.7±0.3 


7 


5401 


2803.531 


Mgll 


0.926±0.001 


1.1±0.3 


8 


5500 


2852.964 


Mgl 


0.928±0.001 


0.6±0.3 






Weighted mean 




0.927±0.001 




9 


3599 


1670.787 


AllI 


1.154±0.002 


1.3±0.4 


10 


3995 


1854.716 


Aim 


1.154±0.002 


0.8±0.3 


11 


5060 


2344.214 


Fell 


1.158±0.001 


0.7±0.3 


12 


5133 


2382.765 


Fell 


1.154±0.001 


0.6±0.3 


13 


5573 


2586.650 


Fell 


1 \ C A \ C\ f\f\ 1 

1.154±0.001 


1.1±0.4 


14 


5603 


2600.173 


Fell 


1.155±0.001 


1.0±0.4 


15 


6024 


2796.352 


Mgll 


1.154±0.001 


0.7±0.3 


16 


6038 


2803.531 


Mgll 


1.154±0.001 


1.2±0.4 






Weighted mean 




1.155±0.001 




i / 




1 707 
1 jZO. /U / 


olll 


1 .4 / oztU.UUZ 


l.ZztU. J 


1 K 


JO JO 


1 ^sn 11 A 

i JJU. / /'-r 


V.1 V 


i .H- / JZLu.UUZ 


Z. JZtlu.Z 


1 Q 


41 JZ 


1 A7n 7Q7 
iO /U. /O / 


AItt 


1.4/ JztU.UUZ 


1. IzhU. J 


20 


4472 


1808 013 


Sin 


1.474±0.002 


0.7±0.3 


21 


5014 


2026.136 


Znll/Crll blend 


1.475±0.001 


0.6±0.2* 


22 


5105 


2062.664 


Znll/Mgl blend 


1.475±0.001 


1.1±0.3 


23 


5800 


2344.214 


Fell 


1.474±0.001 


2.2±0.2 


24 


5880 


2374.461 


Fell 


1.476±0.001 


1.9±0.2 


25 


5894 


2382.765 


Fell 


1.474±0.001 


1.7±0.2 


26 


6384 


2576.107 


Mnll 


1.478±0.001 


0.5±0.2 


27 


6401 


2586.650 


Fell 


1.474±0.001 


0.8±0.3 


28 


6423 


2593.731 


Mnll 


1.476±0.001 


2.0±0.2 


29 


6435 


2600.173 


Fell 


1.475±0.001 


1.9±0.2 


30 


6454 


2605.697 


Mnll 


1.477±0.001 


1.1±0.3 


31 


6919 


2796.352 


Mgll 


1.474±0.001 


2.8±0.2 


32 


6937 


2803.531 


Mgll 


1.474±0.001 


3.4±0.3 


33 


7069 


2852.964 


Mgl 


1.478±0.001 


1.0±0.3 






Weighted mean 




1.475±0.001 





* These lines might be blended with each other 



To be consistent with magnitude measurements appeared 
on the GCN circulars archive^ calibration was done using the 
V BVRI magnitudes, as measured by Henden (2001c), of field 
star 'A' indicated by Stanek et al. (2001a). However we used 
other field stars to check the stability of this calibration: we 
found it to be accurate to within 5%. Any significant short-term 
variation of star 'A' can be ruled out from our data. Due to the 
U — B color dependence of the DOLoReS CCD response in the 
U band, a ^0.2 mag color-term correction was subtracted from 

^ GCN Circulars are available at: 

http: / /gcn.gsfc.nasa. gov/ gen/ gen 3^r chive . html 



the U magnitudes of the OT. We remark that the photometry 
errors quoted throughout the rest of the paper are only statisti- 
cal and do not account for any possible (most likely very small) 
zero point offset. 

We also retrieved and reduced the two y-band images ob- 
tained by BilHngs (2001) on 2001 Feb. 22.502 and 22.547 UT. 
From these, we obtain V = 18.80±0.09 and 19.11±0.11, re- 
spectively, assuming the calibration by Henden (2001c) as de- 
scribed above. 

Standard procedures were used to reduce the NIR data: a 
sky estimate for each image was computed from the clipped 
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Fig. 1. (Left panel) Field of GRB010222 as imaged by the TNG on Feb 23.219 UT in the R band. The OT is indicated by the 
tickmarks. (Right panel) The same field as it appears on the DSS-II sky survey. No source is detected at the OT position. Both 
fields have a size of about 4'x4'; North is at top and East is to the left. We also indicated the star 'A' by Stanek et al. (2001a), 
which we used for the calibration 




Fig. 2. Smoothed i^s-band image acquired at TIRGO on 2001 
Feb 24.1 17 UT and centered at the NIR afterglow position. The 
image size is about 2.'5x2.'5. North is at top, east is to the left. 
The circle indicating the position of the NIR transient has a 
radius of 15" 



median of the nearby images. The telescope dithering was 
measured from the offsets of field objects in each image and 
the images were averaged together using inter-pixel shifts. 
Magnitudes were measured inside circular apertures of diam- 



eter 10" and corrected to total magnitudes. The photometric 
calibration was performed using NIR standard stars from the 
hst by Hunt et al. (1998). 

Given the poor seeing, the objects within 10" from the OT 
location could contribute to the measured NIR flux of the tran- 
sient. It is not possible to accurately measure this effect, but we 
estimate that it is at most a small fraction of the total. Indeed, 
the peak of the emission is in good agreement with the expected 
position of the OT, and no evidence of flux is seen at the po- 
sition of an object detected in the optical at about 5" north of 
the OT. Furthermore, the PSF profiles of the OT and of the ob- 
ject located 13" south-east of the OT itself are well distincted, 
suggesting that the latter does not contaminate the OT flux (see 
Fig. 2). 

We rescaled the TIRGO Kg magnitude to the standard K 
band. This was done in two steps: we first considered that, 
given the almost identical characteristics of the Kg and K' fil- 
ters (Wainscoat & Cowie 1992) we could assume K' = Kg 
within the large errors associated with our Kg measurement. 
Second, using the relation K' ~ K = 0.2 x {H — K) by 
Wainscoat & Cowie (1992) and assuming the power law spec- 
tral shape described in Section 3.3 to estimate the iJ-band mag- 
nitude of the OT, we obtained K' — K ~ 0.15. Thus, K = 
17. 35 ±0.3 at the time of the TIRGO observation. 

We then evaluated the Galactic absorption in the optical 
and NIR bands along the direction of GRBO 10222 using the 
Galactic dust infrared maps by Schlegel et al. (1998); from 
these data we obtained a color excess E{B ~V)= 0.022. By 
applying the law by Rieke & Lebofsky (1985), this colorexcess 
corresponds to Ay = 0.07; then, using the relation by Cardelli 
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et al. (1989), we derived Au = 0.12, Ab = 0.09, Ar = 0.05, Aj 
= 0.04, ^j=0.02, Ak = 0.01. 

Spectra, after correction for flat-field and bias, were back- 
ground subtracted and optimally extracted (Horne 1986) using 
Helium-Argon lamps were used for wavelength cali- 
bration; spectra were then flux-calibrated by using the spec- 
trophotometric standard Feige 34 (Massey et al. 1988) and fi- 
nally averaged together The correctness of the wavelength and 
flux caUbrations was checked against the position of night sky 
lines and the photometric data collected around the epoch in 
which the spectra were acquired, respectively. The typical er- 
ror was 0.5 A for the wavelength calibration and 5% for the 
flux calibration. 



3. Results 

3. 1. Photometry 

We clearly detect the transient in all optical and NIR filters; Fig. 
1 shows the TNG i?-band image taken on 2001 Feb. 23.219 
UT along with the same portion of sky as it appears on the 
DSS-II survey. The OT, indicated by the tickmarks, is clearly 
evident on the TNG image, while no object is present in the 
DSS-II image at that position. The association of this transient 
with GRBO 10222 is confirmed by the fading behaviour of this 
source (see optical magnitudes reported in Table 1). The light 
curves in UBVRI bands are reported in Fig. 3, where our data 
are complemented with those published by other authors (GCN 
circulars archive; Sagar et al. 2001; Staneket al. 2001c; Cowsik 
& Bhargavi 2001). Note that we rescaled the BVRI magni- 
tudes appeared in the GCN circulars to the zero point measured 
by Henden (2001c); we did not plot in Fig. 3 the measurements 
in GCN circulars which could not be tied to this photometric 
zero point (for instance, unfiltered or independently calibrated 
magnitudes). No correction has been applied for Galactic ex- 
tinction, which is anyway small (see previous Section), to the 
data in Table 1 and Fig. 3; nor has any host galaxy continuum 
emission been subtracted, this being unknown at present. 

We find that a single power law decay of the form F oc 
(i — io) where to is the GRB start time, does not provide a 
satisfactory fit of the light curves. Indeed, considering a decay 
index a ^ 0.9 as reported by Price et al. (2001b) and Fynbo 
et al. (2001b), we obtain a fit with xt ^ 65 over 68 degrees of 
freedom (dof). With this fit, the measurements acquired start- 
ing about 1 day after the GRB fall systematically below the 
expected decay. This is clearly evident in the B-, V- and R- 
band light curves, the best sampled since the first hours after 
the GRB. 

We thus hypothesize that the i?-band light curve can be 
modeled with a smoothly broken power law as the one used 
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by Stanek et al. (1999) to fit the optical light curves of the 
GRB990510 afterglow. The model we apply is; 



Fu{t) 



(1) 



where tt is the time of the light curve break, ai is the asymp- 
totic decay index for t <^ ti, and a2 is the asymptotic decay 
index for t ^ tb (times are expressed since to). Note that for 
the temporal decay indices we use inverted notations with re- 
spect to those of Stanek et al. (1999). We obtain a best fit of 
the i?-band data with ai = 0.23±0.01, as = 1.57±0.01 and 4 
= 0.54±0.01 days after the GRB start time; however, this fit is 
not satisfactory, as it has xt ~ 14 (66 dof) and systematically 
underestimates the OT flux for t > 5 days after the GRB. 

We next try the generalization adopted by Beuermann et al. 
(1999) of a smoothly broken power law model: 



F^it) ^ {F,it)-^ + F^ityn- 



(2) 



where Fi{t) — kit^"^ (with i = 1,2), and in which the pa- 
rameter s indicates the smoothness of the change from a decay 
index to the other; for very large values of s, this function as- 
sumes a broken power law shape. This model fits much better 
{X^ = 1.6, 65 dof) the R points than model of Eq. (1). An F- 
test shows that the model by Beuermann et al. (1999) leads to a 
tiny chance probability of improvement with respect to model 
ofEq. (1). 

The formal fit parameters are: ai = 0.54±0.08, a2 = 
1.31 ±0.02, tb = 0.44±0.05 days after the GRB start time, and 
'smoothness parameter' s = 3.1 ±1.0. The value of this parame- 
ter indicates that the slope change occurred quite fastly around 
12 hours after to. However, the fit parameters ai and tb are very 
sensitive to the value of s, and we have noted that similarly ac- 
ceptable fits ixf, ranging from 1.6 to 1.7) are obtained when 
s varies from ^3 to ^10. CoiTespondingly, ai varies between 
^-^0.5 and ^^0.8, and tb between ^0.4 and ^0.7 days. The upper 
bounds of the intervals for these parameters are more consistent 
with the results of other authors (Sagar et al. 2001 ; Stanek et al. 
2001c). Therefore, we will assume for ai and tb their average 
values within the above ranges, and the associated dispersions 
as uncertainties: ai = 0.65±0.15 and tb = 0.55±0.15 days. The 
parameter 0:2 is instead very weakly sensitive to the value of s. 

The xt values associated with the fits are larger than 1 and 
therefore formally not completely satisfactory. By applying a 
systematic error of 1 .5 % to all R data points, in addition to the 
statistical and calibration uncertainties, the xt becomes '-^l. 
Such a small additional uncertainty can be ascribed to non- 
homogeneity of i?-band data set (acquired at different tele- 
scopes and analyzed using different methods), although small 
irregularities in the circumburst medium may produce in the 
afterglow light curve small time scale fluctuations of similar 
amplitude (e.g. GRB000301C, Masetti et al. 2000b; see also 
Wang & Loeb 2000). 

The U B,V, I and J light curves, albeit less well sampled, 
are consistent with the model describing the _R-band data (see 
Fig. 3); therefore, no appreciable color evolution is present in 
the optical afterglow. On the contrary, the iiT-band light curve 
is consistent with being constant between the two epochs of 
observation. 
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Fig. 3. UBVRIJK light curves of GRB010222 afterglow, based on the data presented in this paper and in the literature (see 
text). Different symbol styles indicate different bands. Filled symbols represent data presented in this work, while open symbols 
refer to measurements published by other authors. We have consistently referred all optical magnitudes to the calibration zero 
point of Henden (2001c). No Galactic extinction correction, nor host galaxy flux subtraction have been applied. The ORB start 
time corresponds to 2001 February 22.30799 UT. Overplotted to the R-hsaA data (soUd curve) is the best-fit double power law 
empirical model described in Sect. 3.1. Fit parameters are: ai = 0.65±0.15, a2 = 1.31 ±0.02, if, = 0.55±0.15 days. The same 
model curve, with unchanged model parameters, except for the flux normahzation at the break time, accounts well for the data 
in the other bands (dashed curves), except K 



From the fits, we obtained the following average color 
indices (not corrected for Galactic absorption): U — B = 
-0.55±0.05, B-V = 0.43±0.05, V-R = 0.36±0.05, R-I 
= 0.55 ±0.05 and / - J = 0.72±0.15. These figures place the 
OT of GRB010222 in flie locus populated by GRB optical af- 
terglows in the optical color-color diagrams, as iUustrated by 
Simon et al. (2001). 



3.2. Spectroscopy 

We detect several absorption features in the optical spectrum of 
the GRB010222 afterglow (Fig. 4). As earlier reported by other 
authors (Garnavich et al. 2001; Jha et al. 2001a,b; Bloom et al. 
2001), these correspond to three different metal absorption sys- 
tems located at different redshifts. Our line fitting, performed 
with the SPLOT task within IRAF, assumes a Gaussian profile 



for the absorption lines. A conservative error of ±3 A, com- 
parable with the pixel size, is associated with each line wave- 
length measurement. 

On the averaged TNG spectrum, presented in Fig. 4, we 
mark the positions of all lines we identified. Each number cor- 
responds to a fine in Table 2, where line wavelengths, identifi- 
cations and redshifts are reported. From our line identifications 
we measure the three absorption systems at z = 0.927±0.001, 
z = 1.155±0.001 and z = 1.475±0.001, the highest one thus 
being the lower limit for the redshift of this GRB. These results 
are consistent with those of Garnavich et al. (2001), Jha et al. 
(2001a,b) and Bloom et al. (2001). Our spectral resolution is 
not sufficient to confirm the presence of a fine structure in the 
farthest absorbing system as reported by Castro et al. (2001). 

In Table 2 we also list the equivalent widths (EWs) of the 
identified fines computed in the absorber rest frame (i.e. by di- 
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Fig. 4. TNG+DOLoReS spectrum of the GRB010222 afterglow in the 3500-7700 A range, smoothed with a Gaussian filter with 
(7=3 A (i.e. comparable with the spectral dispersion) and corrected for Galactic absorption assuming E{B — V) = 0.022. Numbers 
mark the positions of identified lines as listed in Table 2. Telluric absorption lines at 6870 and 7600 A are also marked with the 
symbol ® . The dashed curve indicates the position of the spectral continuum we used to measure the EW of fines 



viding the measured value by the factor I+2:). The errors on the 
EWs are computed by assuming different spectral continuum 
levels roughly corresponding to la variation of the continuum 
itself in proximity of each line. 

It can be noted from Table 2 that the metal lines of the ab- 
sorption system at the highest z are very strong. In particular, 
the detection of Mgl suggests that the optical emission of the 
GRB afterglow pierced through a dense medium, most likely 
that of a host galaxy at z= 1 .475 . If we consider the study of 2: < 
1.65 Mgll systems in QSO spectra by Rao & Turnshek (2000) 
we find that, according to these authors, Mgii systems with 
rest frame EW of MgiiA2796 > 0.6 A are primarily associated 
with Damped Lyman-a Absorption (DLA) systems. DLAs are 
characterized by being mostly Hi interstellar medium clouds 
of high redshift galaxies. Indeed, about 50% of the Rao & 
Turnshek's (2000) sample of Mgll absorbers with MgllA2796 
and FeiiA2600 EWs larger than 0.5 A have Hi column density 
Nm > 2 X 10^° atoms cm"^ and basically aU have A^hi > 10^^ 
atoms cm~^. 

Moreover, they found that almost all DLAs with A^hi > 
2 X 10^0 atoms cm'^ have MgiA2852 EW larger than 0.7 
A. Thus, given the EWs of Mgl, Fell and Mgll reported in 
Table 2, we have good indications that the 2=1.475 absorber, 
if not associated with the GRB itself, most likely resembles 
the properties of a DLA system with column density A^hi > 
2 X 10^° atoms cm~^. The ratio of the EWs of the Mgii dou- 
blet is Wrest ( A2796)/ Wrest (A2803) = 0.8±0.1; this indicates 
that the absorption is saturated and, given the low resolution 
of the spectrum, the colutim density is undetermined. Rao & 
Turnshek (2000) also show that there is no precise relation- 
ship between the Hi column density and the rest frame EW 
of MgiiA2796. Thus, although we have good indications that 
the z=lA75 system bears substantial absorption, the Hi column 
density is basically unknown. Applying the same approach to 
the other two absorption systems lying along the GRB line of 



sight (although in the case of the absorber at z=1.155 no Mgl 
is detected), we suggest that these systems have a Hi column 
density in excess of 10^^ atoms cm^^. All this points to the 
possibility of the presence of additional extragalactic absorp- 
tion along the GRBO 10222 fine of sight. 

Finally, we tested the hypothesis (Sari et al. 1998) that the 
optical spectrum can be described with a power law in the 
form F;/ oc . li we &t the TNG spectrum with this law, 
we obtain /3 = 1.16±0.05 (with xt ~ 1), consistent with the 
fit results from the multiwavelength spectra shown in the next 
Subsection, but somewhat steeper than the value reported by 
Jha et al. (2001b). 

3.3. Multiwavelength spectra 

In order to study the temporal evolution of the multiwavelength 
spectrum of the GRBO 10222 afterglow, we have selected 5 
epochs from ^0.2 to ^5 days after the high-energy event, at 
which best spectral sampling is available. We corrected the 
afterglow optical/NlR magnitudes for the Galactic absorption 
and converted them into fluxes using the tables by Fukugita et 
al. (1995) for the optical and by Bersanelfi et al. (1991) for the 
NIR. If strictly simultaneous observations are not available, we 
reduced the fluxes to the reference epoch by interpolating the 
empirical model fitted to the i?-band light curve (see Sect. 3.1), 
and normalizing to the appropriate band. A 7% systematic er- 
ror is added quadratically to the measurements to account for 
interpolation and magnitude-to-flux conversion uncertainties. 
The results are plotted in Fig. 5. 

The optical spectral distributions covering B, V, R and, 
when available, U and / filters are well fit at all epochs by 
a single power law of average spectral index (3 = 1.1 ±0.1 (the 
x1 is typically ^1). A constant spectral shape is consistent with 
the achromatic evolution of the OT decay. 
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Fig. 5. Optical/NIR broad-band spectra of the GRB010222 af- 
terglow at five different epochs: 0.19, 0.97, 1.80, 2.95, and 4.96 
days after the GRB. The data are corrected for Galactic ab- 
sorption. A 7% systematic error was added quadratically to 
the measurement errors in order to account for interpolation 
and magnitude-to-flux conversion uncertainties. For the sake 
of comparison, we also plotted the TNG optical spectrum re- 
ported in Fig. 4 (acquired at epoch = 1.04 days after the 
GRB) 

A single power law fit to the UBVRIJK data at the third 
epoch gives (3 = 1.3±0.1 with = 2.3 (5 dof), shghtly steeper 
than the optical slope. Although this may suggest some absorp- 
tion within the host galaxy in addition to the Galactic one (e.g., 
Lee et al. 2001), we note that the difference is barely signifi- 
cant, considering the uncertainty on the joint optical/NIR spec- 
tral fit. Therefore we conclude that the optical spectra at all 
epochs and the optical/NIR spectrum at 1 .8 days after the GRB 
do not deviate significantly from a single power law of slope 

On the contrary, the X-band point at the second epoch (Af 
= 0.97 days after the GRB) strongly deviates from this trend. 
If we attempt to fit all UBVRIJK points at this epoch with a 
single power law we obtain an unacceptable fit. 

4. Discussion 

The optical light curves of the GRB010222 afterglow are de- 
scribed by a continuously steepening trend, instead of a single 
power law, a behaviour seen in about 1/3 of the well mon- 
itored optical afterglows (GRB980519, Jaunsen et al. 2001; 
GRB990123, Castro-Tirado et al. 1999; GRB990510, Stanek 
et al. 1999; GRB990705, Masetti et al. 2000a; GRB991208, 
Castro-Tirado et al. 2001; GRB991216, Halpem et al. 2000; 
GRB000301C, Masetti et al. 2000b, Jensen et al. 2001; 
GRB000926, Fynbo et al. 2001a). In the present case, this can- 



Fig. 6. Multiwavelength X-ray/optical/NIR SEDs of the 
GRBO 10222 afterglow at epochs 0.97 and 1.8 days after the 
GRB. The X-ray data are from in 't Zand et al. (2001). Optical 
and NIR data were corrected both for Galactic extinction and 
for intrinsic host galaxy absorption using the extinction curve 
by Calzetti (1997) and a local color excess E{B - V) = 0.35. 
The resulting optical/NIR spectrum has a slope (3 ~ 0.55 at 
both epochs (the K-bwd point of the first epoch was excluded 
from the fit because highly deviant). This plot indicates that Vc 
lies between the optical and X-ray ranges, at ~3.5 x 10^^ Hz 



not be ascribed to the synchrotron cooling frequency Uc cross- 
ing the optical band in a spherical regime (Sari et al. 1998), 
because the expected variation of decay index would be in this 
case much smaller (Aa = 0.25) than observed (Aa ~ 0.75). 
Moreover, although the limited optical band hampers a neat 
distinction, the transition should be wavelength dependent, 
while the observed steepening is achromatic. 

Therefore we have considered a decelerating and sideways 
expanding jet scenario (Sari et al. 1999). In this case, the index 
a2 should equal the index of the electron energy distribution, 
thus p= 1.3. From this follows that the optical decay index be- 
fore the colhmation break time must be a\ = 0.23, in case Uc 
is above the optical band at those epochs. This is rather dif- 
ferent from our value of ai, 0.65±0.15. Therefore, within the 
present picture, at the observing epochs preceding the colh- 
mation break time, should have already crossed the optical 
band, so that the predicted ai would be ^0.5, more similar to 
its observed value. The slope of the optical spectrum should 
then be pl2 ~ 0.7, significantly flatter than both our measured 
value /3 = 1 . 1 and that {J3 = 0.9) reported by Jha et al. (2001b). 

Local absorption, on top of the Galactic extinction, may 
account for the difference (see below), but would lead to an 
intrinsic optical spectral slope rather different than found in X- 
rays (Px — 1; in 't Zand et al. 2001), which is inconsistent with 
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the fact that both bands should be above Uc. In addition, the in- 
dex of the electron energy distribution derived in the sideways 
expanding jet scenario, p= 1.3, is very flat (in general, p values 
larger than 2 are found; see Frontera et al. 2000) and may pose 
energetic problems (although one can assume that the electron 
power law has a cut-off at some energy). Emission from the 
stiU undetected host galaxy of GRBO 10222 may flatten flie late 
epoch light curve of the afterglow, but the most recent -R-band 
photometry (our Fig. 3; Stanek et al. 2001c; Sagar et al. 2001) 
suggests that the total magnitude of the host caimot be brighter 
than i? ~ 27 and that its contribution at the latest monitor- 
ing epochs (^40 days after the GRB) is not significant, so that 
subtracting it from the optical measurements does not alter the 
fitted value of a2, and therefore of p. 

Given the problems posed by the application of the above 
model to our optical/NIR afterglow data, we have considered 
two alternatives. 

First we have compared our results with a model of an 
external shock expanding in a pre-ejected wind, which deter- 
mines a dependence for the medium density (Panaitescu & 
Kumar 2000). In the optical/NIR bands, it is predicted that, at 
the early epochs (up to 0.1 days after explosion in the source 
rest frame), the temporal decay is significantly flatter (ai = 
0.25) than found by us, while at later epochs our measured a2 
can be reproduced with p = 2.1, which is acceptable. 

Then, we have made the hypothesis that the temporal break 
observed in the optical light curves may be generated by the 
transition from relativistic to non-relativistic conditions in the 
shock (Dai & Lu 1999). By applying the prescriptions of Dai & 
Lu (1999) to our measured decay indices we obtain compatible 
conditions on p only if the Uc is above the optical band during 
our monitoring. From ai = 0.65±0.15 and a2 = 1.31 ±0.02 we 
derive p = 1.9±0.2 and p = 2.27±0.01, respectively, which we 
consider marginally consistent with p ~ 2.1-2.2. This appears 
to be an acceptable value for the electron distribution shape. 
The predicted spectrum, for z^opt < I'c, is /3 ~ 0.6, much flatter 
than observed. We have made the hypothesis that intrinsic ab- 
sorption at the source produces the observed steeper spectrum 
and we have searched for a de-extinction curve to correct for 
tiiis. 

As noted by Jha et al. (2001b), the optical spectrum lacks, 
at rest frame z = 1.475, the typical 2175 A dust absorption wide 
feature seen in the reddening curve of the Milky Way (Pei 1992, 
Cardelli et al. 1989) and of star forming galaxies (e.g. Calzetti 
1997). Thus, intrinsic absorption must not be large. We have 
considered both an extinction curve appropriate for starburst 
galaxies (Calzetti 1997) and the SMC curve (Pei 1992), and 
applied them to our data. With the latter curve we cannot repro- 
duce the predicted spectral slope in the optical/NIR wavelength 
range for any value of E{B — V) (while in the optical range 
only, acceptable fits are obtained, see also Lee et al. 2001). The 
starburst curve provides instead an acceptable correction and 
reproduces power-law spectra (except for the deviating if-band 
point at epoch 0.97 days) with index P ~ 0.5 for a moderate 
color excess, E{B — F) ~ 0.35 (see Fig. 6). The corrected 
optical/NIR spectral index is much flatter than that found in 
X-rays, which suggests that Vc lies between the optical and X- 



ray bands. This is strengthened by the consistency of the X-ray 
spectral index with p/2. 

An empirical value of Vc — 3.5 x 10^'' Hz (supposed to 
be nearly constant during the non-relativistic expansion, Dai 
& Lu 1999) is determined through the extrapolation of the op- 
tical and X-ray spectra (Fig. 6). The hypothesis of extinction 
within a heavily star forming host galaxy has the advantage of 
removing the mismatch between the optical and X-ray spec- 
tra normalizations (otherwise attributable to other causes, like 
inverse Compton emission in the X-rays, as suggested by in 
't Zand et al. 2001), and had been applied in a similar context 
also to GRB971214 (Dal Fiume et al. 2000). 

The ,7-band light curve behaviour is similar to that of the 
optical ones, indicating that the injection break frequency is 
below this band at both epochs of NIR observations. However, 
the behaviour of the /C-band light curve deviates significantly 
from that observed at shorter wavelengths. Comparison of the 
broad-band spectra at 0.97 and 1.8 days after GRB could sug- 
gest that the injection break may be located between the K and 
J bands at the former epoch, and has moved below the K band 
at the latter. Although this would be consistent with the evolu- 
tion time scale of ly^a (oc t~^^^), the higher fluxes detected at 
sub-mm wavelengths at preceding and successive epochs, indi- 
cate that the injection break should be located longward of the 
NIR band during our monitoring (Fich et al. 2001; Kulkarni et 
al. 2001). Therefore, we have no straightforward explanation 
for the deviation of the ii'-band point from a single optical-to- 
NIR power-law at 0.97 days. 

Atz= 1.475, the total isotropic energy emitted in the 2-700 
keV range by tiiis GRB is 7.7 x 10^^ gj-g (in 't Zand et al. 2001). 
With this total energy, and assuming a circumburst medium 
density of n = 10^ — 10^ cm"-^, the epoch of the transition 
from ultrarelativistic to non-relativistic conditions is about 1 
day (Dai & Lu 1999). Therefore, the most viable interpretation 
of our data appears to be a shock undergoing a substantial de- 
celeration and transition to non-relativistic conditions at about 
0.5 days after the GRB (see also Dai & Cheng 2001). The con- 
sequent requirement of a thick ambient medium (denser than 
10^ cm~^) is supported by the observation of the absorbing 
system at z = 1.475, likely coinciding with the redshift of the 
source. This density would be typical for starburst galaxies 
(Calzetti 1997) or DLAs (Wolfe et al. 1986; Rao & Turnshek 
2000), and is similar to that found for a number of GRB hosts 
(GRB990510, Vreeswijk et al. 2001; GRB000301C, Jensen et 
al. 2001; GRB000926, Fynbo et al. 2001a). Therefore, as for 
many previous OTs, there is the suggestion that GRB010222 
occurred in a site of substantial star formation. 

Our observations and conclusions point to the importance 
of early optical/NIR monitoring of GRB afterglows to establish 
their behaviour and to get insight into the GRB progenitor. This 
wiU be guaranteed by the current (HETE-2) or soon-to-fly GRB 
missions {AGILE and Swift). 
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